Comment on "Orbital-selective Mott transitions in two-band Hubbard models 
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A. Liebsch 

Institut fur Festkorperforschung, Forschungszentrum Julich, 52425 Jiilich, Germany 

A recent paper by BKimer et al. cond- mat/0609758 again criticizes earlier QMC/DMFT results 
by Liebsch [Phys. Rev. B 70, 165103 (2004)]. This criticism is shown to be unfounded. 



In Ref. 1 Bliimer et al. continue to criticize earlier 
QMC/DMFT calculations by Liebsch @ for the non- 
isotropic two-band Hubbard model. Now it is claimed: 
"We quantify numerical errors in earlier QMC data which 
had obscured the second transition" and: "The second 
transition is lost in the noise of earlier data [2] with er- 
rors exceeding 100 % at both transitions" . 

We point out that Ref. 1 once again does not provide 
any comparisons of self-energies or spectral distributions 
[3], nor does it refer to recent work [4,5] which confirms 
the results of Ref. 2. A direct comparison demonstrates, 
as we show here, that both QMC calculations are in good 
agreement and that the above claims are unfounded. 

Fig. 1 shows the comparison of QMC self-energies cal- 
culated in Ref. 2 with those of Ref. 1. Fig. 2 shows analo- 
gous results obtained via exact diagonalization (ED) and 
numerical renormalization group (NRG) @. 
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Evidently, all calculations give the same trend: When 
the narrow band becomes insulating, the self-energy of 
the wide band no longer exhibits ~ iuj n behavior at low 
frequencies, as would be characteristic of a Fermi-liquid. 
Instead, it shows progressive bad-metallic behavior, ap- 
proaching a finite value in the iuj n — > limit. This value 
grows with increasing U, until it diverges near 2.7 eV. 
Ref. 2 states: "^(i^n) becomes inversely proportional 
to oj n at 2.7 eV, i.e., a gap opens up." 

Precisely this behavior is seen in the quasi-particlc 
spectra derived in Refs. 1 and 2 (see Fig. 3). Despite 
the differences caused by different maximum entropy fit- 
ting parameters, the low-frequency region is in perfect 
agreement. Both spectra show that, when the narrow 
band becomes insulating, the wide band reveals a pseu- 
dogap which gets progressively deeper with increasing U, 
until this band becomes fully insulating near 2.7 eV. 
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FIG. 1: (a) QMC subband self-energies Ej(io! n ) for different 
U, from Fig. 10 of Ref. 2. Solid red curves: narrow band; 
dashed blue curves: wide band, (b) QMC self-energy of wide 
band, adapted from Fig. 2 of Ref. 1 (KBD). The narrow band 
is insulating in this range. 



FIG. 2: (a) ED subband self-energies for the same parameters 
as in Fig. 1(a). Solid red curves: narrow band; dashed blue 
curves: wide band, (b) Comparison of ED and NRG self- 
energies of wide band. The narrow band is insulating in this 
range. (Both results from Ref. 5). 
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FIG. 3: Quasiparticle spectra from Ref. 2 at T — 31 meV (left 
panels) and from Ref. 1 (KBD) at T = 25 meV (right pan- 
els). Solid red curves: narrow band; dashed blue curves: wide 
band; dotted curves: bare densities of states. (See Ref. 3.) 

The self-energies and spectral distributions indicate 
that there is good agreement between the QMC results 
of Refs. 1 and 2. Both describe identical physics and are 
consistent with the ED and NRG results. There is no ev- 
idence of any sort of disparity. Thus, the claims in Ref. 1 
have no basis. 

Fig. 4 compares Zj(C/) = 1/[1— lvaSi(wo)/wo\- Again, 
there is good agreement, consistent with the results in 
Figs. 1-3. The main difference is that the Coulomb en- 
ergies at which subbands become insulating are slightly 
lower in Ref. 1 than in Ref. 2 (see also Fig. 3). 

Although the Zi(U) are derived from the self-energies 
in Fig. 1, Ref. 1 claims that AZi(U) (obtained by sub- 
tracting results from Refs. 1 and 2) reveals a qualitative 
difference: "second transition lost in noise" , etc. 

A proper analysis of AZi(U) should, of course, include 
(i) the different error margins resulting from QMC sta- 
tistical uncertainties, number of sweeps and time slices, 
and (ii) the different Coulomb energies at which sub- 
bands become insulating, for instance, as a result of a 
different U mesh, different number of iterations and crit- 




FIG. 4: Comparison of Zi(U) derived within QMC for T = 
31 meV, Ising exchange with J — (7/4 Q- Z\\ narrow band; 
Z2: wide band. Solid and open dots (blue): results of Ref. 2 
(L); crosses (red): results of Ref. 1 (KBD). 

ical slowing down. These issues are particularly impor- 
tant when Ei(iw„) becomes singular and Zi(U) becomes 
small. Since all of this is ignored in Ref. 1, it is no sur- 
prise that the agreement seen in Fig. 4 can be turned, at 
specific points, into ficticious disagreement of AZi(U) of 
arbitrary magnitude. 

Evidently the criticism in Ref. 1 is based on a fun- 
damentally inadequate analysis of AZi(U). Moreover, 
Ref. 1 does not provide the reader with direct compar- 
isons of self-energies or spectral functions, such as given 
here in Figs. 1,3 or in Ref. 3, which demonstrate good 
agreement. Finally, Ref. 1 ignores that the QMC results 
of Ref. 2 were fully confirmed by ED and NRG calcula- 
tions [4,5]. 

We conclude that the QMC/DMFT results of Ref. 2 
are correct: The non-isotropic two-band Hubbard model 
with Ising exchange exhibits a first-order Mott transition 
near U = 2.1 eV when the narrow band becomes insu- 
lating, with characteristic hysteresis behavior, and there 
is no sign of first-order behavior when the wide band be- 
comes insulating near 2.7 eV. 

As also shown by the ED/DMFT calculations in Ref. 4, 
to obtain sequential first-order Mott transitions, it is 
essential to go beyond Ising exchange and include full 
Hund's coupling. 
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